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Abstract 


The Multi-angle Imaging SpectroRadiometer (MISR) Research Aerosol algorithm makes it 
possible to study individual aerosol plumes in considerable detail. From the MISR data for two 
optically thick, near-source plumes from the spring 2010 eruption of the Eyjajjallaj okull volcano, 
we map aerosol optical depth (AOD) gradients and changing aerosol particle types with this 
algorithm; several days downwind, we identify the occurrence of volcanic ash particles and 
retrieve AOD, demonstrating the extent and the limits of ash detection and mapping capability 
with the multi-angle, multi-spectral imaging data. Retrieved volcanic plume AOD and particle 
microphysical properties are distinct from background values near-source, as well as for over- 
water cases several days downwind. The results also provide some indication that as they evolve, 
plume particles brighten, and average particle size decreases. Such detailed mapping offers 
context for suborbital plume observations having much more limited sampling. The MISR 
Standard aerosol product identified similar trends in plume properties as the Research algorithm, 
though with much smaller differences compared to background, and it does not resolve plume 
structure. Better optical analogs of non-spherical volcanic ash, and coincident suborbital data to 
validate the satellite retrieval results, are the factors most important for further advancing the 
remote sensing of volcanic ash plumes from space. 

1. Introduction 

Satellite observations can play a key role in constraining aerosol transport models used to 
diagnose the environmental impacts of volcanic eruptions [e.g., Stahl et al., 2011; Heinold et al., 
2011], Emitted-aerosol microphysical properties are among the most important volcanic plume 
characteristics for air traffic safety, and are also significant indicators of eruption style and 
intensity. The ability to distinguish non-spherical volcanic ash from spherical water and sulfate 
particles near-source [e.g., Scollo et al., 2012], to identify ash concentrations downwind, and to 
constrain particle size, are key contributions multi-angle, multi-spectral remote sensing can 
make, at least in principle, toward characterizing volcanic eruptions. 

This paper explores the ability to retrieve and to map, with data from the NASA Earth Observing 
System’s Multi-angle Imaging SpectroRadiometer (MISR) instrument, particle properties for 
both near-source and downwind volcanic plumes. MISR flies aboard the Terra satellite, in a sun- 
synchronous, polar orbit that crosses the equator on the descending node at about 10:30 AM 
local time. The instrument measures upwelling short-wave radiance from Earth in four spectral 
bands centered at 446, 558, 672, and 866 nm, at each of nine view angles spread out in the 
forward and aft directions along the flight path, at 70.5°, 60.0°, 45.6°, 26.1°, and nadir [Diner et 
al., 1998], Over a period of seven minutes, as the spacecraft flies overhead, a 380-km-wide 
swath of Earth is successively viewed by each of MISR’s nine cameras. As a result, the 
instrument samples a very large range of scattering angles - between about 60° and 160° at mid 
latitudes, providing constraints on particle size, shape, and single-scattering albedo (SSA), for 
particles between about 0.1 and 2.5 pm in diameter. These views also capture air-mass factors 
ranging from one to three, offering sensitivity to optically thin aerosol layers, and allowing 
aerosol retrieval algorithms to distinguish surface from atmospheric contributions to the top-of- 
atmosphere (TOA) radiance, even over relatively bright desert surfaces. 

Unlike the continuous coverage provided by instruments in geosynchronous orbit such as the 
European Space Agency’s Spinning Enhanced Visible and Infrared Imager (SEVIRI), or the 
once-daily global coverage of broad-swath polar-orbiting instruments such as the NASA Earth 



Observing System’s MODerate-resolution Imaging Spectroradiometer (MODIS) instruments, 
MISR observes a given location once every eight days at the equator, increasing to about once 
every two days near the poles. As such, MISR data are complementary to the more frequent data 
sources, uniquely offering periodic, regional-scale maps of aerosol type that can discriminate 
spherical smoke and pollution particles from non-spherical particles, such as desert dust [Chen et 
al., 2008; Kalashnikova and Kahn, 2006] and thin cirrus [Pierce et al., 2010], and providing 
two-dimensional maps of near-source plume height from multi-angle stereo [Moroney et al., 
2002; Muller et al., 2002; Kahn et al., 2007; Scollo et al., 2010], 

The MISR Version 22 Standard Level 2 aerosol retrieval product reports aerosol optical depth 
(AOD) and aerosol type at 17.6 km resolution, based on analyzing TOA radiances from 16x16 
1.1 km-pixel regions, in the context of 74 mixtures of up to three aerosol components [Diner et 
al., 2006; Martonchik et al., 2009; Kahn et al., 2010], For eruptions of Mt. Etna between 2000 
and 2008, Scollo et al. [2012] demonstrated that the MISR V22 Standard aerosol product was 
able to detect even low concentrations of volcanic ash in the atmosphere, and on about ten 
occasions for which there were ground-based validation data, MISR reliably distinguished 
sulfate- and/or water-dominated from ash-dominated plumes. This distinction is indicative of 
eruption strength, and is otherwise unmonitored for most volcanoes around the world. 

The current paper focuses on MISR observations of Eyjajjallajdkull volcanic plumes between 14 
April and 16 May 2010. Table 1 summarizes the MISR observations of the volcanic plume, 
including those that imaged the volcano itself or its immediate surroundings during the study 
period, and several downwind cases when MISR observed locations where the ash plume was 
expected, based on aerosol transport modeling. 

Section 2 reviews our approach to characterizing aerosol amount and type with the MISR 
Aerosol Research Retrieval. In Section 3, we present Research Retrieval results for two near- 
source plumes (Figures la and lb), demonstrating the degree to which aerosol type can be 
constrained from MISR observations of relatively thick aerosol layers, and showing qualitatively 
how plume particle properties evolve within the first few hundred kilometers of the source. 
Section 4 looks at two downwind locations, where the volcanic aerosol plume is more diffuse 
(Figures lc and Id). Here we assess the limits of MISR Version 22 Standard Level 2 and 
Research algorithm retrieval sensitivities to AOD and aerosol type, with more general 
implications for detecting and mapping volcanic aerosols far from their sources. Finally, a 
summary and conclusions are given in Section 5. 


2. Approach to MISR Volcanic Plume Aerosol Type Characterization 

During the month-long study period that covers MISR observations of the spring 2010 
Eyjajjallajdkull eruptions, MISR imaged the volcanic plume near-source 13 times (Table 1), 
eight of which included the volcano itself. Many of these are heavily cloud-covered, in part due 
to ambient meteorological conditions, and in part because volcanic ash itself can nucleate ice 
particles [e.g., Seifert et al., 201 1], We analyze two of the most cloud-free near-source cases in 
detail using with the MISR Research Aerosol Retrieval algorithm [Kahn et al., 2001], and 
compare with MISR V22 Standard algorithm [Martonchik et al., 2002; 2009] results. Whereas 
the Standard algorithm processes the entire MISR data stream automatically, the Research 
algorithm provides considerably more flexibility, at the cost of longer computation time per 
retrieval, as well as significant hands-on preparation required for each individual case. We have 
also expanded some of the data analysis capabilities of the Research algorithm since earlier 



publicaions, and present them for the fist time in this study. Here are key attributes of the 
Research algorithm runs: 

• Retrievals were carried out on equivalent reflectance values averaged, channel-by-channel, 
over retrieval regions either 3 x 3 or 5 x 5 pixels in size. As the near-source retrievals are over 
ocean, angular views within about 30° of the specular glint angle were eliminated to avoid glint 
contamination (40° over regions away from the optically thick plumes). A standard ocean 
surface model was adopted [e.g., Martonchik et al., 1998], with near-surface wind speed 
constrained by NCEP and MERRA reanalysis values [Kalnay et al., 1996 and Rienecker et al., 
2011, respectively], 

• Retrievals were performed over many locations for each event, covering areas defined 
manually, with the aim of characterizing plume properties and avoiding condensate cloud where 
possible. The size of the individual retrieval regions was determined based on scene variability, 
to minimize heterogeneity within retrieval regions while maximizing coverage to the extent 
possible. (The Standard V22 algorithm performs retrievals on a fixed, global grid of 16 x 16 1.1 
km pixel regions.) 

• Over optically thick, elevated plumes, the multi-angle views for each event were manually co- 
registered at the approximate plume elevation (derived from the MISR stereo imagery; Figure 2), 
to minimize camera-to-camera aliasing. 

• Where possible, retrievals were also performed in cloud-free regions adjacent to the plumes, to 
determine background AOD and aerosol type. 

• The retrieval approach is based on selecting aerosol optical models and corresponding AOD 
that produce acceptable matches to the MISR-observed top-of-atmosphere equivalent 
reflectances. Thirty-two aerosol component optical analogs were considered in this study, in four 
aerosol type categories (Table 2): (1) five spherical non-absorbing particles of different sizes, (2) 
20 spherical absorbing particles, representing five sizes at each of four SSA values, all with 
spectrally varying (“steep”) SSA, (3) three non-spherical medium dust grain optical models of 
varying SSA (i.e., modeled as having 1%, 4%, and 10% hematite), plus one coarse dust spheroid 
model, and (4) three non-spherical cirrus optical analogs of different sizes. (The Standard V22 
algorithm considers eight components [Table 2 of Kahn et al., 2010].) The incremental 
differences among the spherical components considered are based on studies that established the 
limits of aerosol type retrieval sensitivity for MISR-like data (e.g., Kahn et al., 2001; Chen et al., 
2008). For practical reasons, the non-spherical “dust grain” and “cirrus” component optical 
models used here are based solely on previous work [Kalashnikova et al., 2005 and Baum et al., 
2005, respectively], 

• Mixing groups were formed, comprised of up to four aerosol components, one each from the 
spherical non-absorbing, spherical absorbing, non-spherical medium dust grains or non-spherical 
coarse-dust spheroids, and non-spherical cirrus categories. In total, 1200 mixing groups were 
included in the comparison space, representing all possible combinations of four-component 
groupings. For each mixing group, retrieval runs tested mixtures of the four components in all 
possible proportions of mid- visible AOD in 10% increments, a total of 286 mixtures for each 
mixing group (and a total of 343,200 mixtures) at each AOD value. (The Standard V22 
algorithm considers a total of 74 specific mixtures of up to three components [Table 3 of Kahn et 
al., 2010].) All tests were performed for total mid- visible AOD values over ranges selected 



based on conditions, in increments of either 0.02 or 0.05, depending on range of AOD to be 
covered. 


• Three normalized tests are used to assess the degree to which modeled values match the 
observed equivalent reflectances, based on comparisons between the measurements and 
simulated values for different choices of AOD and aerosol mixture: (1) absolute reflectances, (2) 
angle-by-angle relative reflectances for each wavelength, and (3) spectral relative reflectances at 
each angle, as described in detail elsewhere [Kahn et al., 2001], For each AOD and mixture in 
the modeled comparison space, the three X tests are performed; the largest value among the 
three tests, representing the largest model-measurement discrepancy compared to uncertainty, is 
designated ma x3- 

• The measurement uncertainty used in calculating X parameters is nominally taken as 5% of the 
observed reflectance, but we vary the acceptance criteria to account for the ambiguity in these 
values. Specifically, we examine results for each of seven different acceptance criteria, three 
absolute criteria: p ^ maX 3 < 3.0, < 2 .0, and < 1.0, and four relative criteria: if p m j n is the minimum 
value of p ^ maX 3 for all mixtures in a given retrieval region, we accept all mixtures in the region 
having ^f maX 3 within p mm + 0.5, /j> nim + 0.25, p mm + 0.1, and the mixture within each retrieval 
region having the smallest ;^ m ax 3 overall. We also examine the channel-by-channel measured 
and model reflectances for individual cases (what we call “signal plots”), to resolve any issues in 
the statistics. Conclusions are drawn based upon the systematic attributes of the AOD and 
aerosol type solution spaces, as more stringent absolute and relative acceptance criteria are 
applied. 

• Given the considerable number of cases involved in exploring the aerosol type parameter space, 
the approach we take for summarizing the results statistically is important. To provide context 
for the numerical values, Figure 3 illustrates the spatial distribution of Research Retrieval results 
for (a) AOD, (b) Angstrom exponent (ANG), (c) fraction AOD non-spherical, (d) single- 
scattering albedo (SSA), (e) fraction AOD “grains,” (f) fraction AOD “cirrus,” (g) fraction AOD 
spherical non-absorbing particles, and (h) fraction AOD spherical absorbing particles, for the 
near-source plume on 07 May 2010, both within the main and secondary plumes, and in a control 
region to the west, superposed on true-color MISR images of the region. These results will be 
discussed subsequently, in Section 3 below. Figure 4 provides several statistical summaries of 
the retrieval results for the nine 5x5 pixel retrieval regions in Patch 2, which is located along 
the main plume, as indicated in Figure 3c. The p ^ maX 3 criterion applied for the data shown is 
given in the upper right of Figure 4, and the total number of mixtures that pass this criterion, 
summed over all retrieval regions in the patch, and all mixing groups included in the summary, is 
reported as the “Count” in the upper left. 

Sensitivity studies indicate that confidence in MISR aerosol-type retrieval results increases when 
mid- visible AOD exceeds about 0.15 or 0.2, and that confidence in identifying individual aerosol 
components increases when the component contributes at least about 20% to the total AOD 
[Kahn et al., 2001; 2010], We use these criteria as general guidelines in constructing and 
interpreting the statistical summary plots, to help extract robust conclusions from the mass of 
retrieval results. The left side of Figure 4 provides the overall-patch and retrieval-region-specific 
Research Retrieval results, taking account of all 1200 mixing groups, and organized by aerosol 
type. 



The uppermost pie chart on the left in Colu mn 1 shows the fractional distribution of broad 
aerosol type categories, aggregated with equal weight over all the retrieval regions contained in 
the patch; the legend below also reports the numerical fraction of each aerosol type reflected in 
the chart. To the right of the aggregated chart are the corresponding charts for each individual 
retrieval region, organized with roughly the geographical spatial distribution, as it appears in 
Figure 3. This makes it possible to identify aerosol-type gradients within the patch, and to assess 
variability. Below these plots on the left side of Figure 4 are similar pie charts - this time 
showing the fractional AOD contributions of the specific aerosol components in those mixtures 
that pass the acceptance criterion. A listing of the components (see Table 2) is given at the 
bottom, along with their fractional contributions to the total AOD, which were first assessed for 
each retrieval region, and then aggregated with equal weight over the entire patch; the values for 
the three components having the highest fractions are shown in bold. 

The remaining columns of Figure 4 report the AOD values and specific, AOD-weighted particle 
property retrieval results, aggregated over all mixtures from all 1200 mixing groups meeting the 
X max3 acceptance criterion. In the lowest three plots, component particle properties are shown, in 
this case including all components contributing to the total AOD, as indicated by the labeling in 
the right margin. (We also look at aggregations that include only those components contributing 
at least, e.g., 20% to the total AOD (AOD/ ;m ., > 0.2), to take advantage of the greater aerosol type 
sensitivity at higher AOD.) Again, overall patch and individual retrieval region results are given 
to the left and right, respectively. In descending order: AOD and ANG, and then the AOD 
fractions assigned to spherical vs. non-spherical, and to different ranges of r e and SSA. Not 
shown are similar figures for the other six x?max 3 acceptance criteria, and for particle properties 
where AOD y ract > 0.2 instead of AOD y ract > 0.0. 

In addition to the near-source events presented in Section 3, we studied two cases where satellite 
imagery and HYSPLIT [ Draxler and Rolph, 2003] trajectory modeling indicated likely plume 
occurrence far downwind of the volcano, but within the MISR field-of-view. Whereas the near- 
source analysis allows us to look at variability and changing particle properties within the plume, 
the downwind case analyses focus on the limits of MISR ability to detect and map volcanic 
plumes. 


3. Near-Source Volcanic Plume Particle Properties from MISR 

Volcanic aerosol component and mixture emissions can vary on time scales as short as hours or 
less, and subsequently, plume aerosol properties can evolve as the particles age and the plume 
dissipates. Often, especially near-source, volcanic plumes contain a distinct super-micron 
(coarse) mode, dominated by silicate ash that can vary in composition, size, shape, and 
absorption properties, and a sub-micron (fine) mode, dominated by sulfuric acid and/or 
ammonium sulfate, that is generally spherical and weakly or non-absorbing; the proportions of 
coarse to fine-mode can also vary in space and time, and such behavior is observed for the 2010 
Eyjajjallajdkull eruptions [e.g., Schumann et al, 201 1; Ansmann et al., 2011], MISR observations 
are snapshots of plumes, giving an instantaneous view of particles that are progressively older 
downwind, and offering some information about particle microphysical properties, generally: 
three-to-five size bins, two-to-four bins in SSA, and spherical vs. non-spherical shape [Kahn et 
al., 2001; 2010], A main strength of these data is that they provide extensive spatial coverage, 
offering loose constraints on plume evolution, and context for more detailed suborbital 
measurements. 



3.1. 07 May 2010 Near-source Plume Properties 


At 12:39 UTC on 07 May 2010, MISR imaged the Eyjajjallajdkull plume, starting at the volcano 
itself and reaching -550 km downwind (Figures la). The observations capture a range of 
conditions from a relatively strong eruption (Table 1) - fresh emissions near the source to 
transported aerosol about 6-10 hours old, based on forward-trajectory analysis [Draxler and 
Rolph, 2003], MISR stereo heights report a plume vertical extent from approximately 4 km to 6 
km above the ocean surface near-source, and beginning about 400 km downwind, it descends to 
2-3 km along with the ambient air parcels (Figure 2a). MISR stereo-derived, near-source plume 
vertical-extent characterization for the Eyjajjallajdkull volcano is presented in more detail, and 
assessed by comparison with other published plume height observations, by Garay et al. [2012], 

The MISR Standard Level 2 Aerosol product provides few results directly over the plume core, 
in part because the algorithm does not co-register the images to the height of the plume itself, 
and in part due to high AOD and scene variability over the 17.6 km scale at which retrievals are 
performed by the algorithm. The adjacent retrieval regions, however, show distinctly plume-like 
characteristics compared to background values. The retrieved mid- visible AOD is - 1.1, about 
twice that of the surroundings, and the particle ANG - 0.3, SSA - 0.95, and fraction AOD non- 
spherical - 0.7 - 0.8. By comparison, values retrieved for the background are ANG > 0.75, SSA 
> 0.97, and fraction non-spherical <~ 0.4, i.e., smaller, brighter, and more spherical particles, as 
might be expected. 

The plume broadened horizontally downwind, and became more horizontally uniform. Figure 3 
shows the MISR Research Retrieval results geographically distributed over the plume, and Table 
3 summarizes the results for four patches along the main plume (outlined in Figure 3c), plus a 
low-altitude secondary plume of remobilized ash very near the volcano itself (Patch S), and a 
control region over relatively clear ocean just west of the main plume (Patch C). 

In Figure 3, the plume aerosol stands out as distinct from the background, having much higher 
AOD, and containing much larger (lower ANG), far larger fraction non-spherical, and generally 
darker (lower mid-visible SSA) particles than the control region to the southwest, as might be 
expected. Overall, mid-visible AOD in the plume core exceeds 4 from near-source to about 400 
km downwind, ANG is below -0.25, and SSA is around 0.93-0.95, based on the MISR Research 
retrieval results. In comparison, the MISR Standard retrievals reflect similar plume particle 
property trends, but the differences relative to background are much smaller, and the plume 
structure is not resolved. 

Peak AOD decreases systematically beyond about 250 km downwind of the source (Table 3), 
and retrieved aerosol properties along the plume edges increasingly converge toward background 
values. In those parts of the plume beyond 400 km that were imaged by MISR, as the plume 
begins to descend, AOD drops considerably. 

From about 250 km downwind to 500 km along the plume, retrieved grain and cirrus analogs 
vary, but in sum, these non-spherical species consistently contribute about 65% to the total 
retrieved AOD (Table 3). Similarly, in the plume core, the AOD contribution attributed to small, 
spherical, non-absorbing particles remains fairly constant, at -10-15% throughout the plume 
(Figure 3g), as does the AOD fraction of spherical absorbing particles, at -20-25%. However, 
the SSA assigned to the spherical absorbing particles increases from -0.80 near-source to -0.86 
550 km downwind, and the retrieved median ANG increases a small amount, from about 0.12 to 



0.20 (Table 3), which might indicate gradual brightening and decreasing particle size, 
respectively, as the plume ages. 

As discussed subsequently, mixtures of the spherical absorbing and the non-spherical 
components probably represent the optical properties of volcanic ash in these retrievals. Note 
that although the retrieved “cirrus” optical component for the patches over most of the main 
plume is identified with the smallest-size option (10pm, Table 2), for the very near-source 
Secondary Plume, about 25% of the cirrus is attributed to the larger (40 pm) component (Table 
3). Similarly, the Secondary Plume retrieval allocates a larger AOD fraction to spherical 
absorbing particles than any of the other retrieval regions, suggesting these are darker particles 
and/or contain a larger fraction of volcanic ash. 

Actual cirrus cloud (meteorological cirrus) is ubiquitous in the area, possibly nucleated in places 
by the plume itself, though also associated with what appears to be a frontal boundary just to the 
east of the field-of-view, that can be seen more clearly in the corresponding, larger-swath 
MODIS/Terra image (not shown). The cirrus is easily observed in the background image of 
Figure 3d, which shows the MISR 70° forward image rather than the nadir view. Taking a closer 
look at the plume structure, Figures If, h, and j show cirrus dominating the NE edge from PI to 
~15 km upwind of P2, and again from P2 to the MISR swath edge near P3. From the multi- 
angle views, the cirrus is elevated with respect to the plume between PI and P2, which precludes 
successful retrievals at some locations within the plume. Figure 3e indicates a dominance of 
grains over the same stretch of the plume core; apparent variations in grain fraction could 
represent fluctuations in volcanic emission content during the eruption, though difficulty in co- 
registering the multi-angle views due to multiple layers at different elevations might also be 
involved. Looking downwind to P3 and P4, about 400 to 550 km from the source, the grain 
contribution remains high (>30-45%) over most of the plume extent. At the southernmost edge 
of P4, there is indication that cloud contamination is again important. 

Figure 4 focuses on the retrieved aerosol properties for P2, 250 km downwind from the source. 
(See Section 2 above for a detailed explanation of the features of Figure 4.) The nine retrieval 
regions in this patch were chosen to represent the plume core. The retrieval results for less 
stringent acceptance criteria are similar, and converge systematically toward the values shown in 
Figure 4, as described in Section 2. (The plots for other acceptance criteria are not shown.) 

The left side of Figure 4 makes it possible to parse the retrieved properties shown in Figure 3 
into component particles at P2. There is notable consistency in the component particle types 
derived for the individual retrieval regions in the patch, though the proportions vary. Overall, 
more than 60% of the plume AOD is attributed primarily to non-spherical grains and 10 pm 
cirrus optical analogs, about a fifth to small-medium, spherical, absorbing species having SSA 
around 0.8, and about 14% to small-medium, non-absorbing spherical particles having effective 
radius (r e ) mostly 0.12 pm. Focusing on the non-spherical components, nearly half the total 
AOD is attributed to medium dust grains modeled as having 1% or 4% hematite (weakly or 
moderately absorbing, Table 2), and about 15% to 10 pm cirrus. 

Referring to the 1200-mixing-group comparison space used for this study, spherical absorbing or 
non-absorbing particles larger than 0.57 pm did not match well, and neither did ellipsoids (Table 
2). Very small (;y=0.06 pm) particles were not selected in early retrieval runs, and were therefore 
eliminated from the final comparison space used for the figures. However, the identification of 
specific aerosol components based on the MISR retrievals in this situation is complicated by a 
lack of good optical analogs for the volcanic ash particles [e.g., Schumann et al., 2011], likely 



cirrus contamination, scene heterogeneity, and difficulty co-registering some parts of the scene. 
If the atmospheric columns at the plume core were composed of spherical, non-absorbing sulfate 
particles plus volcanic ash, typical of volcanic plumes and confirmed by in situ observations of 
Eyjajjallajdkull plume particles at other times and places, it is likely that ash particles are 
represented optically in the retrievals by combinations of weakly to moderately absorbing non- 
spherical particles combined with small-medium spherical absorbing particles. The water/sulfate 
droplets could be represented by the small spherical non-absorbing particle component and 
possibly some fraction of the spherical absorbing particles. Further interpretation, based on 
comparison between these retrieval results and suborbital measurements of the 2010 
Eyjajjallajdkull plume particles, is given in Section 3.3 below. 


3.2. 19 April 2010 Near-source Plume Properties 

Based on surface observations, the 19 April eruption was less energetic than that the 07 May 
event (Table 1), so lower AOD would be expected, and possibly also lower ash/spherical particle 
ratio and smaller mean particle size. The plume is shown in Figure lb; note the sharp edge along 
the western side of the ash plume for about the first 100 km from the source. As this plume is 
less elevated and less textured than the 07 May case, the MISR Standard Level 2 aerosol product 
obtained retrievals over the plume itself on 19 April. Mid- visible AOD is about 1.0, about twice 
the background value. The retrieved plume particles are distinctly larger than background, with 
ANG -0.1 compared to >0.75, and there is a somewhat larger fraction non-spherical, ~0.8 
compared to -0.6; no distinct pattern appears in the MISR Standard Retrieval product SSA. 
Figure 5 provides an overview of the Research Retrieval results, and Table 4 gives retrieval 
summaries for six patches progressively downwind. 

For over-water retrievals, we usually exclude cameras that view within 40° of the specular or 
“sun-glint” direction to avoid un-modeled brightness when interpreting the TOA reflectances. 
The differences in available cameras based on this criterion extend along vertical strips through 
the entire image. This plume is not as optically thick as the 07 May case, so glint from the ocean 
surface can make a larger contribution to the TOA reflectances; at the same time, this plume is 
extensive in the east-west direction, so the number of cameras removed by the glint constraint 
varies from two along the western side of the plume to four at the eastern edge. When the 
retrievals were performed with varying numbers of cameras by the Research algorithm, the range 
of observed scattering angles varied between about 45° and 90° from east to west, so the 
constraints on aerosol type varied, though not uniformly along the camera difference lines, and 
the corresponding retrieved AOD also changed discontinuously in places, by as much as 60% in 
the higher AOD, near-source region. So the retrievals shown in Figure 5 consistently use five 
cameras over the entire plume, to eliminate artifacts that having different numbers of cameras 
can produce in this situation. 

Retrieved particle microphysical properties are more variable, and the plume appears to evolve 
much more rapidly that in the 07 May case; the AOD diminishes downwind, with peak retrieved 
values exceeding 3.0 only within about 20 km of the source, decreasing systematically to <2.0 
for patches beyond 50 km (Table 4). The variations in the retrieved snapshot of plume properties 
might be due in part to plume dissipation and the aging of plume particles with observed distance 
from the source, and changes in the erupted material over the 9-12 hours represented by the 
volcanic effluent shown in Figure 5, based on trajectory modeling [Draxler and Rolph, 2003], 
Other possible factors contributing to the variability include condensate cloud, which is common 
in this scene (e.g., Figure 5d). In addition, MISR stereo heights suggest the plume is 



concentrated between 0.25 and about 1.5 km ASL, yet within about 100 km of the source, some 
1.1 km stereo-height retrieval elements of the plume occur as high as 4 km (Figure 2b). 


Nevertheless, retrieved microphysical properties are distinct from background maritime aerosol, 
and characteristic of volcanic plumes. From beyond 20 km downwind to 275 km, retrieved 
particle properties remain fairly uniform (Table 4); small, spherical non-absorbing particles 
comprise 14±3% of the total AOD based on the retrieval results, and grains contribute about 
47±3%. Unlike the 07 May case, the grains are mostly highly absorbing (10% hematite), 
suggesting that this eruption might have produced generally darker material. Cirrus analogs add 
about another 15±4% to the total AOD within the plume, and in this case are dominated by the 
large (40+100 pm) components. Note also that outside the plume, just west of PI in Figures 5e 
and 5f, the AOD fraction contributed by non-spherical grains and cirrus is near zero, the AOD 
drops below 0.5 (Figure 5a), and the ANG is above 0.75 (Figure 5b), providing a significant size 
and shape contrast to the particles within the volcanic plume. Again the MISR Standard 
retrievals reflect similar but more muted plume particle property tendencies, and much less 
structure. 

Beyond about 325 km downwind, particle properties begin to change, suggesting that the aging 
process might become more important. For Patch P6, 400 km downwind, a larger fraction of the 
retrieved grains is weakly (1% and 4% hematite) rather than highly absorbing, the retrieved 
cirrus analogs are mostly small (10 pm) instead of large, and the SSA for the mix of particles 
increases from -0.88 near the source to -0.91 - a small absolute SSA difference, but more robust 
as a relative indication [e.g., Kahn et al., 2009], Taken together, these results seem to support 
the idea that as the volcanic plume ages, the average aerosol properties tend toward smaller and 
brighter particles. A lack of coincident, in situ validation data limits the confidence with which 
we present these interpretations, but in the next subsection, we glean what we can about the 
particle microphysical properties from available aircraft measurements of the spring 2010 
Eyjajjallajdkull volcano eruptions. 


3.3. Comparisons with Suborbital Particle Property Measurements 

Although no aircraft in situ samples of plume particles were collected coincident with the MISR 
observations, the DLR Falcon F20 obtained two near-source samples from this eruption 
sequence, on 02 May over the North Atlantic within 7-12 hours of injection, and on 17 May over 
the North Sea, 60-84 hours after injection [Schumann et al., 201 1], The analysis by Schumann et 
al. emphasizes the largest airborne ash particles and determination of layer extent and particle 
concentration, due to the relevance of these quantities for aircraft safety. As MISR has greater 
sensitivity to the properties of particles about 0.1 pm and 2.5 pm in diameter [Kahn et al., 2010], 
our comparisons focus on particles in this size range. 

Chemical and morphological analysis of their aircraft samples by Schumann et al. [201 1], from 
both 02 and 17 May, identified a mix of ash and sulfate particles, with particles smaller than 0.5 
pm in diameter dominated by sulfate, and those larger than 0.5 pm dominated by silicate ash, 
though the proportions of ash/sulfate and the mineralogical composition of the ash varied among 
the samples. Nearly all the light-absorbing particles fell in the size range 0.5- 1.0 pm diameter, 
and for the samples acquired on 02 May, this size range corresponds to ash mixed with a 
significant fraction of sulfate. As the imaginary index-of-refraction is uncertain from the 
analysis, hematite-containing ash was assumed to be the absorbing species in all cases, and 



interpretation of airborne optical measurements in terms of particle size distributions was treated 
parametrically. 


For the best-estimate imaginary index-of-refraction, nearly all observed particles had diameters 
below 20 pm, as soon as 10 hours after injection and beyond just a few hundred km downwind 
of the source; the peak diameter of the coarse mode fell between 2.2 and 13.5 pm, and the fine 
mode peaked at about 0.1 pm diameter [Schumann et al., 2011], The fine mode was interpreted 
as nucleating sulfate particles. Best-estimate overall effective diameter ( D e ff= 3V/2A, where V 
and A are the particulate volume and cross-sectional area per unit volume of atmosphere, 
respectively) spanned the range of 0.2 to 3 pm, and did not show any trend with plume age 
among the aircraft measurements. Regarding particle absorption, the best-estimate effective 
SSA at 532 nm, derived from T-matrix calculations of ellipsoidal particles initialized with the 
aggregate of aircraft-constrained effective size distributions and an assumed imaginary refractive 
index of 0.001, yields values around 0.95 for particles smaller than 0.5 pm, 0.9 for those between 
0.5 and 1 pm, descending to about 0.85 for particles having A# that exceeds about 2 pm. Note 
that larger A# in this dataset correlates with higher mass concentration, which in turn 
corresponds to observations made nearer the volcanic source [Schumann et al., 2011], 

Placing the MISR Research Retrieval results in the context of these in situ observations, the 
effective size of the spherical, non-absorbing fine-mode particles are in good agreement. 
Comparison of the coarse-mode, and more generally, the ash particles, is more difficult for a 
number of reasons: (1) there are no suborbital measurements coincident with MISR observations, 
and the coarse mode optical properties in particular are highly variable in the volcanic plume in 
both space and time, (2) MISR sensitivity to particle properties (though not AOD) is weighted 
toward particles having diameter <~ 2.5 pm whereas the aircraft analysis was weighted toward 
larger particles due to a focus on aircraft safety, (3) particle index-of-refraction is poorly 
constrained in both data sets; this significantly affects interpretation of particle size distribution 
in the aircraft data and also contributes to MISR retrieval uncertainty, and (4) particle shape was 
not measured by the aircraft instruments, and interpretation of the aircraft data used ellipsoids 
optical models for the ash, whereas the MISR retrievals rejected ellipsoid optical models in favor 
of more angular grain and “cirrus” models. 

Given the aggregate of uncertainties, and in particular, the apparent lack of a good optical model 
for the volcanic ash, it is likely that the ash particles are represented optically in the MISR 
retrievals in most cases by combinations of (1) medium, weakly to moderately absorbing non- 
spherical grains, (2) non-spherical ~10 pm cirrus, and (3) small-medium, spherical absorbing 
particles, in proportions that vary based on plume evolution and retrieval ambiguities. Two key 
conclusions from the analysis in this section are that MISR is sensitive to the general 
microphysical properties of near-source volcanic ash plumes, and that there is a pressing need for 
additional detailed, direct, suborbital measurements of volcanic ash optical properties to help 
constrain remote sensing retrievals. 


4. MISR Volcanic Plume Particle Detection and Mapping Far Downwind 

The near-source retrieval studies presented in Section 3 above focused on gleaning from the 
MISR data as much information as possible about the spatial distribution of plume physical and 
optical properties. This section treats two far-downwind cases, more than a day subsequent to 
atmospheric injection. It aims primarily at assessing the limits of plume detection, based on 



combinations of AOD, particle shape, and to some extent particle size contrasts with background 
aerosol. 


4.1 .16 April 2010 Plume Detection Near Cabauw 

Figures lc and 6 show MISR images of the Eyjajjallajdkull volcanic plume near Cabauw at 
about 10:45 UTC on 16 April, approximately 24-36 hours downwind of the source [Draxler and 
Rolph, 2003], The cloud-enshrouded coast of the Netherlands is barely visible just east of the 
relatively cloud- free water over which Research Retrievals were performed. The brown ash 
plume is most easily seen in Figure 6d, where the background image is the MISR long-slant-path 
70° forward view. In the MISR data, the plume is generally thinner than the near-source cases, 
but it remained fairly narrow (~ 60 km wide, depending on what AOD is used to define the 
plume edge), and the maximum mid-visible AOD in the plume core is still ~1 .35. The MISR 
Standard aerosol product did not provide results in this region, due to the shallow water mask; 
however, the imagery (Figure lc) and Research retrieval results suggest the top-of-atmosphere 
reflectances are atmosphere-dominated in this case. MISR stereo height retrievals indicate a 
plume vertical distribution extending between 1 and 3 km ASL, with most of the plume 
concentrated between 1.5 and 2.5 km. 

Although AOD decreases toward the plume edges, particle properties remain fairly constant, as 
might be expected for a gradually dissipating feature. Retrievals in the plume core (Patch P in 
Figure 6c) assign about half the AOD to non-spherical, weakly absorbing grains containing 1% 
or 4% hematite. High grain content relative to background is a distinctive feature of the volcanic 
ash plume retrievals for all the cases in this study, though trajectory analysis lends confidence to 
the volcanic plume identification for the far-downwind observations. About 20% of the AOD in 
this case is attributed to small-medium, non-absorbing spherical particles, about 15% to small- 
medium spherical absorbing particles having effective SSA around 0.86, and -10% to large 
spherical absorbing particles having SSA of -0.80. (Figure 6 and Table 5). Very small, spherical 
particles, as well as coarse-mode ellipsoid particles are excluded in these retrievals, and unlike 
the near-source cases, only 5% of the AOD is associated with cirrus analogs. For the aggregate 
of particles, the retrieved ANG falls between about 0.1 and 0.5, and the SSA is between 0.90 and 
0.95. 

Compared with the near-source plumes discussed in Section 3, this aged plume is again 
represented by small-medium, spherical non-absorbing particles that could be optical analogs of 
sulfate or water particles, plus a mixture of weakly absorbing non-spherical and medium-large 
spherical absorbing particles that combine to represent volcanic ash optically. The main 
difference is that the downwind plume has a much smaller fraction of the larger non-spherical 
cirrus analogs, which suggests that the effective size of the non-spherical (ash) component 
decreased as the plume evolved. 

This event was also detected, beginning around 12:00 UTC on 16 April, by the European 
Aerosol Research Lidar Network (EARLINET) station at Cabauw [Ansmann et al., 2011]. The 
site is located about 57 km inland (52.0 N, 4.9 E) from the offshore MISR retrieval region. As 
there is considerably more cloud cover over the land (Figure 6), the AERONET station at the 
Cabauw site [ Holben et al., 1998] reported highest-quality (Level 2) results only beginning four 
hours after the MISR overpass. At that time, the AERONET AOD was descending rapidly to 
below 0.3, from values probably in excess of 0.5 two hours earlier, based on available Level 1.5 
data. The Cabauw lidar identified an ash layer between 2 and 3 km ASL, similar to the nearby 
MISR stereo height determination. 



Although EARLINET particle properties are not reported for this specific event at Cabauw, the 
MISR ANG values between 0.10 and 0.50 are consistent with values of 0.30-0.50 described for 
the transported ash plume at the Leipzig EARLINET site. The MISR-retrieved high fraction of 
super-micron particles (>65% for most areas) and the corresponding low ANG are also 
supported by suborbital data [Ansmann et al., 201 1], Similarly, the mean SSA of the MISR- 
retrieved mixtures appears to match reasonably well the SSA constrained by refractive indices 
from Schumann et al., [2011], though it is lower than the 0.97 value obtained by Hervo et al. 

[201 1] for more aged volcanic plume aerosol over central France on 19 May. Overall, the degree 
to which MISR particle property results coincide with suborbital values appears to substantiate 
the suborbital indications that ash particles >15 pm settled out before reaching continental 
Europe [Ansmann et al., 2011; Flentje et al., 2010; Bukowiecki et al., 2011, though sampling 
biases might preferentially eliminate larger particles]. 


4.2 . 10 May 2010 Plume Detection in the Central North Atlantic 

At 13:13 UTC on 10 May 2010, MISR viewed a plume of transported ash around 5 1 .37 N, 3 1 .06 
W in the central North Atlantic, approximately 1,500 km SSW of the EyjaJjallajdkull volcano. 
The plume was observed here some three days after emission, based on forward and backward- 
trajectory modeling [. Draxler and Rolph, 2003], The MISR Standard aerosol product identified 
higher AOD, lower ANG, and higher fraction AOD non-spherical over the remnant plume than 
in the surrounding region, but the area involved amounts to no more than four 17.6 km retrieval 
regions, so only by comparison with the more detailed and higher-resolution Research Algorithm 
results, discussed below, is it possible to develop confidence associating the results with actual 
plume aerosol properties. MISR stereo analysis could not determine the height of the ash plume 
due to a lack of distinct plume features that must be identified in multiple angular views for the 
technique to work, so the patches in the image were co-registered to 1.25 km, the stereo-derived 
height of the nearby clouds. This assumed height produced consistent retrievals over the plume 
area, lending confidence to the choice. 

The 10 May event is the lowest AOD case in the present study, with peak mid- visible AOD 
-0.85 (Table 5); only the SE comer of the retrieval area, around Patch P in Figure 7c, is clearly 
identifiable in the MISR data as part of the plume. The Control area, Patch C, has AOD in the 
range of 0.15 to 0.2, typical of background values. The Patch P particle properties are distinct 
from the background and similar to upwind volcanic ash plumes, with more than 80% of the 
AOD assigned non-spherical particles - over half to weakly absorbing grains, and nearly a 
quarter retrieved as 10 pm cirrus. The median ANG value is -0.10. In contrast, in the Control 
patch representing background aerosol, no grains are retrieved, and the Angstrom exponent is 
higher, with a median value of -0.37 (i.e., smaller effective particle size), despite having a 
significant AOD fraction (36%) retrieved as large cirrus, which in this case is likely to be actual 
cirrus, based on the meteorological context (Figure Id). 

One difference between this aged plume and the upwind volcanic plumes is that the effective 
SSA for the retrieved spherical absorbing particles is higher here than for the upwind cases 
(Tables 3, 4, and 5). Retrieved SSA for the background particles is lower, but the background 
AOD is so low that SSA retrieval results are not well constrained [e.g., Kahn et al., 2010], SSA 
is expected to increase as plumes age, due to oxidation and coating of ash particles, though 
validation data to confirm the satellite remote-sensing result is lacking for this eruption. 



5. Conclusions 


We have presented a detailed analysis of AOD and particle microphysical properties for two 
near-source and two far downwind MISR observations of the spring 2010 Eyjajjallajdkull 
volcano plume. The MISR retrievals provide maps of aerosol amount and type within the 
volcanic plume, showing plume structure and evolution, at least qualitatively. They complement 
more limited suborbital measurements, as well as more frequent but less detailed broad-swath 
polar-orbiting and geostationary satellite instrument observations, by providing basic plume 
property information where it is otherwise lacking. Case studies such as these are also key steps 
at assessing the limits of MISR retrieval capabilities, and can serve as templates for volcanic 
plume analysis procedures that might subsequently be applied elsewhere around the globe. The 
MISR Research Aerosol algorithm makes it possible to select specific retrieval region locations 
and sizes, to constrain surface reflectance properties based on time-series analysis, to consider 
hundreds of aerosol components and hundreds of thousands of mixtures, and to explore a range 
of acceptance criteria options. It is the main tool used in this study. 

We find that both retrieved volcanic plume AOD and particle microphysical properties are 
distinct from background values near-source, as well as for over- water cases several days 
downwind, for situations where the plume remained relatively concentrated. Cloud 
contamination for this particular series of events, and the general challenge of adequately 
characterizing the surface reflectance over land, precluded clear identification of volcanic ash 
over land sites several days downwind, which we tested at several sites where EARLINET lidar 
detected thin ash layers aloft (Table 1). Especially for the far-downwind cases, trajectory 
analysis played an important role in identifying regions where the volcanic plume was likely to 
reside. 

The MISR Standard Version 22 aerosol product is not optimized for this application due to 
coarse spatial resolution, a highly constrained number of aerosol components especially 
regarding spherical, absorbing particles) and mixtures in the algorithm climatology, and a lack of 
image co-registration at individual plume elevations. Within these limitations, the standard 
product did obtain higher AOD, and distinctly larger, darker, and more non-spherical particles 
compared to background values, all trends generally expected, and reflected in the Research 
algorithm results. However, the MISR Standard retrievals derived much smaller differences 
compared to background than the Research algorithm results, and did not resolve plume 
structure. 

Where detected, the volcanic ash particles were characterized optically as mixtures of non- 
spherical grains and cirrus analogs, plus small-medium spherical absorbing particles, due to a 
lack of good optical models for ash itself. Small-medium, spherical non-absorbing particles, and 
possibly some part of the spherical absorbing component AOD, are probably associated with 
sulfate and water particles typically present in volcanic plumes, and identified in non-coincident, 
suborbital observations of this eruption series. As 866 nm is the longest wavelength available 
with the MISR instrument, retrieval sensitivity to particle microphysical properties (but not the 
AOD) decreases for particles larger than about 2.5 pm in diameter, though there is still some 
ability to distinguish larger cirrus analogs of different sizes [Pierce et al., 2010], In situ 
observations suggest that the largest size for transported ash from spring 2010 the 
Eyjajjallajdkull eruptions was about 15 pm. The MISR Research retrievals precluded very small 
particles (<0.12 pm effective radius), and available ellipsoidal ash optical analogs (Table 2); they 
showed strong sensitivity to the differences between weakly and highly absorbing grains, and 



between 10 pm and 40-100 pm cirrus analogs. The results also offer qualitative indication of 
expected trends in particle properties as the plumes aged: particle brightening and decreased 
average particle size, showing to some degree the spatial and temporal patterns of plume particle 
evolution. 

Unlike the multi-year time series of Mount Etna observations [ Scollo et al., 2012], there was no 
discemable pattern of ash-to- sulfate/water particle AOD ratio changes in the few cases studied 
here. This is due at least in part to a lack of good volcanic ash optical models, which depend 
upon having better particle shape, size, and index-of-refraction constraints than those currently 
available. There was also a lack of coincident, in situ particle property data to validate the MISR 
retrievals quantitatively, highlighting the need to acquire such data in the future. As such, the 
results presented here demonstrate in yet another research area where aerosols affect the Earth 
environment, the essential, complementary nature of satellite and suborbital measurement, and 
aerosol transport modeling [e.g., Kahn, 2011]. 
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Figure Captions 

Figure 1. MISR true color context images, showing the two near-source plumes and the two 
downwind plumes from the Spring 2010 eruptions used in this study, (a) 07 May 2010 plume, 
Orbit 55238, Path 216, Blocks 40-42, 12:39 UTC, nadir view, (b) 19 April 2010 plume, Orbit 
54976, Path 218, Blocks 39-42, 12:51 UTC, 70° forward view, (c) 16 April 2010 plume, Orbit 
54931, Path 197, Block 49, 10:45 UTC, 70° forward view, (d) 10 May 2010 plume, Orbit 55282, 
Path 221, Block 49, 13:13 UTC, 70° forward view. Inserts show the approximate geographic 
location of the imagery. 

Figure 2. MISR stereo height retrieval maps derived using the MISR INteractive explorer 
(MINX, Nelson et al., 2008; http://misr.jpl.nasa.gov/getData/accessData/MisrMinxPlumes). (a) 
07 May 2010 plume; Plume 1 is the main plume and Plume 2 is a secondary plume of 
remobilized ash. (b) 19 April 2010 plume. All heights are above sea level (ASL), and are 
reported at 1 . 1 km horizontal resolution. A wind correction has been applied, and the vertical 
accuracy of the resulting height retrievals is approximately 500 m. 

Figure 3. MISR Research Retrieval results for the 07 May 2010 plume, superposed on MISR 
true-color imagery, (a) AOD, (b) ANG, (c) Fraction AOD non-spherical, (d) SSA, (e) fraction 
AOD “ grains,” (f) fraction AOD “cirrus,” (g) fraction AOD spherical non-absorbing particles, 
and (h) fraction AOD spherical absorbing particles. The volcano is in the upper left of these 
images, north is roughly toward the top, and the nadir image is shown in all but panel (d), where 
the 70° forward view is given, to highlight aerosol detail and thin cloud adjacent to the main 
plume. For scale, the width of the swaths shown is about 380 km. Two plumes appear, the main 
plume at about 5 km elevation, extending to the SE, and a small plume of re-mobilized ash near 
the volcano itself, within 0.5 km of the surface, that is blown to the SW (see Figure 2a). The 
surface wind speed was set to 7.5 m/s, and for the main plume retrievals, the multi-angle images 



were co-registered to the approximate plume elevation. Retrievals were performed on regions 
comprised of 5 x 5 1.1 km pixels over the main plume, the secondary plume, and in a control 
area to the west. The p ^ maX 3 threshold for the retrievals in these images was set to p m i„+ 0. 1 . 
Retrieval results for Patch 2, one of six patches outlined in panel (c), are summarized statistically 
in Figure 4. An overview of Research Retrieval results for all six patches is given in Table 3. 

Figure 4. Summaries of the MISR Research Retrieval results for the nine retrieval regions, each 
5x5 1.1 km pixels in size, located over the main plume on 07 May 2010 in Patch P2 of Figure 
3c, for the AOD-weighted averages of all mixtures in the 1200 mixing groups that meet the 
acceptance criterion ^f maX 3 < p m in + 0.1. (See text, Section 2, for a detailed explanation of the 
features of this figure.) 

Figure 5. Same as Figure 3, but showing results for the 19 April 2010 plume. The volcano is in 
the top center of these images, and north is roughly toward the top. The multi-angle images were 
co-registered to the approximate plume elevation of 3 km, and the surface wind speed was set to 
10 m/s. In this case, five MISR cameras were used consistently over the entire retrieval region, 
to avoid artifacts caused by an increased number of view angles in areas where fewer cameras 
would have been eliminated by the standard glint mask. The p ^ maX 3 threshold for the patches in 
this image was set at p m i„+ 0.1. 

Figure 6. Same as Figure 3, but showing results for the 16 April 2010 plume far downwind near 
Cabauw. North is roughly toward the top. The multi-angle images were co-registered to 2 km, 
the approximate elevation of the ash plume, and the surface wind speed was set to 5 m/s. The 
X max3 threshold for the patches in this image was set at p m ,„+0. 1 . 

Figure 7. Same as Figure 3, but showing results for the 10 May 2010 plume in the central North 
Atlantic. North is roughly toward the top. The elevation of the ash plume could not be 
determined from stereo imaging, so the multi-angle images were co-registered to the nearby 
cloud elevation of 1.25 km; the surface wind speed was set to 2.5 m/s. The ;^ maX 3 threshold for 
the patches in this image was set at p m i n + 0. 1 . 



Table 1. MISR observations of the during the Eyjajjallajdkull plume during the study period. 


Date 

Orbit 

Path 

Block 

Time 

(UTC) 

Notes 

15- April 

54917 

206 

41 

11:38 

Injection height 5-9.5 km; strong westerlies!; 
~1300 km downwind height: 3-5 km**; 
high mass-flow event'*' 

16- April 

54931 

197 

49 

10:46 

Plume ht. 1-3 km near Cabauw** 

18- April 

54961 

211 

41 

12:09 

Less ash; injection ht. 3-5 kmf **; 
high mass-flow event*' 

19-April* 

54976 

218 

39-42 

12:51 

Injection height 0.25-4.5 km**; 
low mass-flow, low ash*; low SC> 2 # 

21 -April 

55005 

216 

40 

12:39 

Weak eruption, ash 1-3 kmf**; 
Low wind, layering! 

24- April* 

55049 

221 

39 

13:10 

Ash layer 1-3 km!; 

~1 .5 km, directed eastward** 

26-April* 

55078 

219 

39 

12:58 

Ash in cloud layer 4.5-5 km** 

3-May* 

55180 

220 

39 

13:04 

Plume darker than in April!; 
Injection height 4±0.5 km** 

5-May* 

55209 

218 

39 

12:51 

Ash > 8km, reached Spain & Morocco, 
strong satellite H 2 S0 4 signal!; 
Injection height 3.5-6 km** 

7-May* 

55238 

216 

40-43 

12:39 

Injection height 2.25-6 km + 
remobilized ash near surface**; 
high mass-flow event*; high SC> 2 # 

10-May 

55282 

221 

46-48 

13:13 

Representative lat-lon (55.16,-28.09); 

~1,000 km SW of source; 

3 days downwind from 5/7 eruption** 

12-May* 

55311 

219 

39 

12:58 

Injection height 3. 5-6 .5 km 

12-May 

55311 

219 

58 

13:04 

Representative lat-lon (40.68,-29.57); 

~2,600 km SSW of source; 

5 days downwind from 5/7 eruption** 

13-May 

55325 

210 

41 

12:03 

After 08 May, plume injection height 



decreased to 6-7 knit; 
Downwind plume ~5 km 5 


Injection height 5-8 km § ; 

16-May* 55369 215 40 12:33 high mass-flow event* 

* Cases where the volcano itself appears in the MISR imagery. 

f Observations by Petersen (2010); plume heights from Keflavik airport radar, about 150 km 
west of EyjaJJallajdkull. 

# SO 2 measurement from the Ozone Monitoring Instrument; available from: 

http://disc.sci.gsfc.nasa.gov/Aura/data-holdings/OMI/omso2_v003.shtml 
J From Table 3 of Schumann et al. (2011). 

§ Injection heights from MISR stereo analysis [ Garay et al., 2011] 

0 Based on HYSPLIT forward and back-trajectories [Draxler and Rolph, 2003] 



Table 2. Aerosol Component Optical Models 


# 

Component 

Name 

n 

(nm) 

r 2 

(urn) 

r c 

(l*m) 

a 

SSA 

(446) 

SSA 

(558) 

SSA 

(672) 

SSA 

(866) 

AOT(446)/ 

AOT(558) 

A ^ 

CD £ 
< < 

AOT(867)/ 

AOT(558) 

g 

(558) 

Particle 

Size/Shape 

Category 


Spherical Non-Absorbing Optical Models* 

1 

sph_nonabsorb 
_0. 1 2 

0.003 

0.75 

0.06 

1.70 

1.00 

1.00 

1.00 

1.00 

1.54 

0.66 

0.35 

0.61 

Small Spherical 

2 

sph_nonabsorb 

_0.26 

0.005 

1.70 

0.12 

1.75 

1.00 

1.00 

1.00 

1.00 

1.18 

0.82 

0.58 

0.72 

Small Spherical 

3 

sph_nonabsorb 

_0.57 

0.008 

3.81 

0.24 

1.80 

1.00 

1.00 

1.00 

1.00 

0.98 

0.99 

0.91 

0.72 

Medium Spherical 

4 

sph_nonabsorb 
_1 .28 

0.013 

8.88 

0.50 

1.85 

1.00 

1.00 

1.00 

1.00 

0.96 

1.04 

1.10 

0.73 

Large Spherical 

5 

sph_nonabsorb 

_2.80 

0.022 

19.83 

1.00 

1.90 

1.00 

1.00 

1.00 

1.00 

0.98 

1.02 

1.05 

0.77 

Large Spherical 


Spherical Absorbing Optical Models* 

6 

sph_abs_0.1 2_ss 
a_green_ 
0.80_steep 

0.003 

0.75 

0.06 

1.70 

0.82 

0.80 

0.77 

0.72 

1.470 

0.695 

0.403 

0.61 

Small Spherical 
very strongly 
absorbing 

7 

sph_abs_0.1 2_ss 
a_green_ 
0.85_steep 

0.003 

0.75 

0.06 

1.70 

0.87 

0.85 

0.83 

0.79 

1.488 

0.686 

0.388 

0.61 

Small Spherical 
strongly 
absorbing 

8 

sph_abs_0.1 2_ss 
a_green_ 
0.90_steep 

0.003 

0.75 

0.06 

1.70 

0.91 

0.90 

0.89 

0.85 

1.507 

0.677 

0.375 

0.61 

Small Spherical 
moderately 
absorbing 

9 

sph_abs_0.1 2_ss 
a_green_ 
0.95_steep 

0.003 

0.75 

0.06 

1.70 

0.96 

0.95 

0.94 

0.92 

1.524 

0.669 

0.361 

0.61 

Small Spherical 
weakly absorbing 

10 

sph_abs_0.26_ss 

a_green_ 

0.80_steep 

0.005 

1.69 

0.12 

1.75 

0.79 

0.80 

0.80 

0.79 

1.165 

0.838 

0.614 

0.75 

Small Spherical 
very strongly 
absorbing 

1 1 

sph_abs_0.26_ss 

a_green_ 

0.85_steep 

0.005 

1.69 

0.12 

1.75 

0.84 

0.85 

0.85 

0.84 

1.171 

0.833 

0.603 

0.74 

Small Spherical 
strongly 
absorbing 

12 

sph_abs_0.26_ss 

a_green_ 

0.90_steep 

0.005 

1.69 

0.12 

1.75 

0.89 

0.90 

0.90 

0.90 

1.176 

0.828 

0.593 

0.73 

Small Spherical 
moderately 
absorbing 

13 

sph_abs_0.26_ss 

a_green_ 

0.95_steep 

0.005 

1.69 

0.12 

1.75 

0.95 

0.95 

0.95 

0.95 

1.181 

0.824 

0.584 

0.73 

Small Spherical 
weakly absorbing 




sph_abs_0.57_ss 

14 a_green_ 0.008 3.805 0.24 

0.80_steep 

sph_abs_0.57_ss 

15 a_green_ 0.008 3.805 0.24 

0.85_steep 

sph_abs_0.57_ss 

16 a_green_ 0.008 3.805 0.24 

0.90_steep 

sph_abs_0.57_ss 

17 a_green_ 0.008 3.805 0.24 

0.95_steep 

sph_abs_1 .28_ss 

18 a_green_ 0.013 8.884 0.5 

0.80_steep 

sph_abs_1 .28_ss 

19 a_green_ 0.013 8.884 0.5 

0.85_steep 

sph_abs_1 .28_ss 

20 a_green_ 0.013 8.884 0.5 

0.90_steep 

sph_abs_1 .28_ss 

21 a_green_ 0.013 8.884 0.5 

0.95_steep 

sph_abs_2.80_ss 

22 a_green_ 0.022 19.826 1 

0.80_steep 

sph_abs_2.80_ss 

23 a_green_ 0.022 19.826 1 

0.85_steep 

sph_abs_2.80_ss 

24 a_green_ 0.022 19.826 1 

0.90_steep 

sph_abs_2.80_ss 

25 a_green_ 0.022 19.826 1 

0.95_steep 


26 dust_grains_h1 0.10 1.00 0.5 

27 dust_grains_h4 0.10 1.00 0.5 


1.80 0.77 0.80 0.82 0.84 0.984 

1.80 0.82 0.85 0.87 0.89 0.982 

1.80 0.88 0.90 0.91 0.93 0.981 

1.80 0.94 0.95 0.96 0.96 0.980 

1.85 0.77 0.80 0.83 0.86 0.962 

1.85 0.82 0.85 0.87 0.90 0.962 

1.85 0.88 0.90 0.92 0.93 0.961 

1.85 0.94 0.95 0.96 0.97 0.961 

1.90 0.77 0.80 0.82 0.85 0.983 

1.90 0.83 0.85 0.87 0.89 0.983 

1.90 0.88 0.90 0.91 0.93 0.982 

1.90 0.94 0.95 0.96 0.97 0.982 

Dust Grains Optical Models* 

1.5 0.92 0.98 0.99 1.00 0.90 

1.5 0.72 0.91 0.98 0.99 0.90 



0.986 

0.915 

0.78 

Medium Spherical 
very strongly 
absorbing 


0.986 

0.914 

0.76 

Medium Spherical 
strongly 
absorbing 


0.987 

0.914 

0.75 

Medium Spherical 
moderately 
absorbing 


0.987 

0.914 

0.74 

Medium Spherical 
weakly absorbing 


1.038 

1.091 

0.78 

Large Spherical 
very strongly 
absorbing 


1.039 

1.093 

0.77 

Large Spherical 
strongly 
absorbing 


1.039 

1.094 

0.76 

Large Spherical 
moderately 
absorbing 


1.040 

1.095 

0.74 

Large Spherical 
weakly absorbing 


1.018 

1.051 

0.83 

Large Spherical 
very strongly 
absorbing 


1.018 

1.051 

0.81 

Large Spherical 
strongly 
absorbing 


1.017 

1.050 

0.80 

Large Spherical 
moderately 
absorbing 


1.017 

1.050 

0.79 

Large Spherical 
weakly absorbing 


1.07 


1.08 


0.71 


Weakly 

absorbing grains 
Moderately 
absorbing grains 


1.07 


1.10 


0.72 



28 

dust_grains_h1 0 

0.10 

1.00 

0.5 

1.5 

0.98 

0.80 

0.94 

0.98 

1.05 

1.09 

1.16 

0.72 

Strongly 
absorbing grains 


Dust Spheroid Optical Model* 

29 

dust_spheroids 

0.10 

6.0 

1.0 

2.0 

0.81 

0.90 

0.97 

0.98 

0.99 

1.02 

1.05 

0.77 

Coarse Dust 
Spheroids 


Cirrus Optical Models* 1 ^ 

30 

Baum_cirrus_ 
De= 1 Oum 

2.0 

9500.0 

5.0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.79 

Cirrus 

31 

Baum_cirrus_ 

De=40um 

2.0 

9500.0 

20.0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.81 

Cirrus 

32 

Baum_cirrus_ 
De= 1 OOum 

2.0 

9500.0 

50.0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.87 

Cirrus 


* Aerosol components are named based on particle shape (spherical, non-spherical grains, spheroids, or cirrus), SSA (non- 
absorbing, weakly, moderately, strongly, or very strongly absorbing) and effective radius (in micrometers). Single 
scattering properties were calculated using a Mie code for the spherical particles; the dust component properties were 
calculated using the Discrete Dipole and T-matrix approaches for medium and coarse modes, respectively [ Kalashnikova et 
al., 2005], Wavelength in nm is specified in parentheses where appropriate, ri and iy are the upper and lower limits of the 
size distribution, r c and a are the characteristic radius and width parameters in the log-normal distribution, and SSA is the 
single-scattering albedo. The asymmetry parameter (g) will generally represent particle scattering phase functions poorly 
for the purpose of calculating MISR multi-angle radiances, and is given here only in MISR green band for reference. The 
vertical distributions of particles are based on plume height; aerosol is concentrated near-surface for control regions. 


t Gamma distribution is used for cirrus particle models, and r c is replaced by r e , the effective radius [Baum et al., 2005]. 




Table 3. Research Retrieval results overview for the 07 May 2010 plume.* See Figure 3c for context. 


Location 
(Fig. 3c label) 

Peak 

AOD 

Median 

ANG 

% AOD 
Grains' 

% AOD 
“Cirrus” f 

% AOD Spherical, 
Non-Absorbing 1 ^ 

% AOD (SSA) Spherical, 
Absorbing 1 ^ 











W 

H 

S 

L 

S 

M 

L 

S 

M 

L 

Secondary Plume (S) 

~ 4.5 

-0.20 

10 

3 

26 

8 

12 

1 

— 

11 (0.83) 

27 (0.82) 

1 (0.80) 

100 km (PI) 

>4.5 

-0.35 

13 

2 

37 

— 

19 

— 

— 

— 

7 (0.80) 

15 (0.80) 

250 km (P2) 

>4.5 

-0.12 

47 

1 

15 

2 

13 

2 

— 

5 (0.82) 

14 (0.83) 

3 (0.82) 

400 km (P3) 

-4.25 

-0.15 

40 

1 

17 

4 

12 

1 

— 

8 (0.86) 

11 (0.83) 

5 (0.81) 

550 km (P4) 

-2.0 

-0.20 

23 

9 

29 

5 

11 

2 

— 

10 (0.86) 

6 (0.86) 

2 (0.82) 

Control Region (C) 

-0.22 

-0.85 

28 

2 

1 

2 

40 

— 

— 

24 (0.90) 

1 (0.86) 

— 


* Results presented here and in the figures are based on the ;f m ax3 < p m ,n +0.1 criterion. 

f The retrieved grains are reported as the percent of the total AOD assigned to weakly or moderately absorbing (W) 
1% and 4% hematite components (Table 2), and the percent of the total AOD assigned to highly absorbing (H) 

10% hematite component. The retrieved percent of the total AOD assigned to smaller (S) 10 pm Cirrus is 
followed by the percent of the total AOD assigned to larger (L) 40+100 pm Cirrus. For the spherical particles, 
percentages of the total AOD and SSA are reported for small (S; <0.30 pm), Medium (M; 0.30<r e <0.70 pm), and 
Large (L, >0.70 pm) components. 

Note that retrieved grains are nearly all the weakly absorbing components, the retrieved “Cirrus” is nearly all 10 pm 
in size, and the spherical non-absorbing particles are effectively all small. The exceptions are (1) the 550 km patch, 
where 9% of the total AOD is assigned to the highly absorbing grains, and (2) the Secondary Plume patch, where 
9% of the total AOD is assigned to 40 pm Cirrus. 




Table 4. Research Retrieval results overview for the 19 April 2010 plume. * See Figure 5c for context. 


Location 
(Fig. 5c label) 

Peak 

AOD 

Median 

ANG 

% AOD 
Grains' 

% AOD 
“Cirrus” + 

% AOD Spherical, 
Non-Absorbing 1 ^ 

% AOD (SSA) Spherical, 
Absorbing 1 ^ 











W 

H 

S 

L 

S 

M 

L 

S 

M 

L 

20 km (PI) 

~ 3.1 

-0.25 

5 

41 

3 

8 

17 

3 

— 

12 (0.88) 

10 (0.86) 

1 (0.81) 

50 km (P2) 

~ 2.0 

-0.20 

4 

40 

6 

12 

12 

2 

— 

15 (0.86) 

7 (0.84) 

2 (0.82) 

75 km (P3) 

~ 1.55 

-0.30 

1 

48 

2 

10 

15 

1 

— 

17 (0.87) 

4 (0.84) 

1 (0.82) 

175 km (P4) 

~ 1.65 

-0.15 

3 

43 

5 

11 

11 

2 

1 

15 (0.85) 

7 (0.83) 

3 (0.83) 

275 km (P5) 

~ 1.2 

-0.25 

4 

45 

6 

10 

15 

2 

1 

13 (0.88) 

3 (0.87) 

2 (0.85) 

400 km (P6) 

~0.8 

-0.20 

17 

30 

14 

7 

11 

1 

1 

12 (0.88) 

4 (0.85) 

2 (0.84) 


* Results presented here and in the figures are based on the )Cmaxi < p m ,n +0.1 criterion, 
f The columns are defined as in Table 3. 




Table 5. Research Retrieval results overview for the downwind plumes. 


Date & 

(Figure, Location) 

Peak 

AOD 

Median 

ANG 

% AOD 
Grains' 

% AOD 
“Cirrus” + 

% AOD Spherical, 
Non-Absorbing 1 ^ 

% AOD (SSA) Spherical, 
Absorbing 1 ^ 











W 

H 

S 

L 

S 

M 

L 

S 

M 

L 














16 April (Fig. 6c, P) 

~ 1.35 

-0.15 

47 

1 

3 

2 

17 

2 

2 

11 (0.86) 

3 (0.87) 

11 (0.80) 














10 May (Fig. 7c, P) 

- 0.85 

—0.10 

58 

— 

23 

— 

7 

2 

1 

5 (0.91) 

2 (0.92) 

1 (0.89) 

10 May (Fig. 7c, C) 

-0.20 

-0.37 

— 

— 

25 

36 

6 

9 

2 

16 (0.83) 

3 (0.83) 

4 (0.84) 


* Results presented here and in the figures are based on the ^f maX 3 < p m ,„ +0. 1 criterion, 
f The columns are defined as in Table 3. 
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